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Gamma-ray emission of relativistic jets as a supercritical process 
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Supercriticality of the same kind as that in a nuclear pile can take place in high-energy astrophysical objects pro- 
ducing a number of impressive effects. For example, it could cause an explosive release of the energy of a cloud of 
ultrarelativistic protons into radiation. More certainly, supercriticality should be responsible for energy dissipation of 
very energetic relativistic fluids such as ultrarelativistic shocks in gamma-ray bursts and jets in active galactic nuclei 
(AGNs). In this case, the photon breeding process operates. It is a kind of the converter mechanism with the high- 
energy photons and e + e _ pairs converting into each other via pair production and inverse Compton scattering. Under 
certain conditions, which should be satisfied in powerful AGNs, the photon breeding mechanism becomes supercritical: 
the high-energy photons breed exponentially until their feedback on the fluid changes its velocity pattern. Then the 
system comes to a self-adjusting near-critical steady state. Monte-Carlo simulations with the detailed treatment of 
particle propagation and interactions demonstrate that a jet with the Lorentz factor T ss 20 can radiate away up to 
a half of its total energy and for Y = 40 the radiation efficiency can be up to 80 per cent. Outer layers of the jet 
decelerate down to a moderate Lorentz factor 2-4, while the spine of the jet has the final Lorentz factor in the range 
10-20 independently on the initial T. Such sharp deceleration under the impact of radiation must cause a number of 
interesting phenomena such as formation of internal shocks and an early generation of turbulence. 
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1. Supercritical processes in 
high-energy astrophysics 

The term "supercriticality" has different meanings. 
Here we use the same meaning as in nuclear physics: 
the condition of exponential runaway behaviour of 
the system (nuclear explosion or nuclear pile) . So far 
two kinds of supercritical phenomena were proposed 
in high-energy astrophysics: an explosive energy re- 
lease of a cloud of ultrarelativistic protons into pho- 
tons ("proton bomb") and dissipation of the bulk 
energy of relativistic fluids into radiation. 

Let us consider a magnetized cloud of relativistic 
protons with typical Lorentz factor j p . These pro- 
tons do not radiate much, because the synchrotron 
radiation and the inverse Compton scattering rates 
for protons are suppressed by factor (m e /m p ) 2 with 
respect to electrons. Suppose N seed photons with 
energy e t ~ 200MeV/7 P interact with protons pro- 
ducing pions decaying into high-energy electrons or 
photons. They, in turn, produce a wide spectrum of 
softer photons, with the number £iV of them on aver- 
age interacting with the high-energy protons again. 
If £ > 1 then the process has a supercritical charac- 



ter: the number of photons rises exponentially until 
the protons cool down, so that the supercriticality 
condition breaks down and the energy release ter- 
minates. If the energy supply to protons continues 
then the system reaches supercriticality and explodes 
again. Such scenario was proposed and studied with 
numerical simulations [1] and analytically [2]. 

The supercriticality condition depends on the 
radiation spectrum emitted by the particles pro- 
duced in photo-meson interaction. The hardest pos- 
sible spectrum is the fast-cooling power-law dNf ds oc 
£ -3/2^ rp-^g so ft es t spectrum is that of the saturated 
pair cascade dN/de oc e~ 2 [3], which takes place in a 
dense soft photon field. 

The cross-section of photo-meson production 
multiplied by the fraction of energy transfered to 
pions, x, is (see [4]) a(p r ))x ~ 0.7 x 10~ 28 cm 2 in 
the photon energy range ~ 250 MeV-2 GeV in the 
proton rest frame. The number of produced photons 
which can interact again is N ~ 1/7 P in the first 
case and N ~ I/7 2 in the second case. Then condi- 
tion £ > 1 translates to pe > 10 25 erg cm -2 in the 
fast cooling case and pe > 10 25 /7 P erg cm -2 in the 
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cascade case, where pe is the column energy density 
of protons across the cloud. 

Where one can find such a cloud of ultrarela- 
tivistic protons? A possible site for this phenomenon, 
gamma-ray bursts (GRBs), was proposed by [5]. In- 
deed, all scenarios of GRBs involve ultrarelativistic 
shocks, where the space between the shock front and 
the contact discontinuity is filled with ultrarelativis- 
tic protons. 

Another possible class of supercritical phenom- 
ena can take place in relativistic flows. A first exam- 
ple is the supercritical pair loading in AGN jets near 
the accretion disc's hot corona which emits hard X- 
rays [6]. If the jet Lorentz factor at such distance is 
T ~ 10, then cool pairs moving with the jet can up- 
scatter X-rays to MeV energies and these photons 
in their turn, interacting with the same X-rays, can 
produce new pairs. 

A more powerful and universal way of dissipa- 
tion of relativistic fluids into radiation, is associated 
with the converter mechanism suggested in general 
form by Derishev et al. [7] and in a more specific form 
of a runaway (supercritical) electromagnetic cascade 
in ultrarelativistic shocks by Stern [8] . We later per- 
formed a numerical study of the converter mecha- 
nism, in its electromagnetic version, operating in rel- 
ativistic jets of blazars [9, 10] and found that it has a 
supercritical character: the high-energy photons and 
relativistic e+e~ pairs breed exponentially. 

2. Converter mechanism and photon 
breeding 

The converter mechanism is an alternative to the 
Fermi acceleration for relativistic flows in a dense 
radiation field. In the Fermi mechanism a charged 
particle gains energy being scattered many times 
between media moving with respect to each other. 
In the converter mechanism, a quantum of energy 
moves between media in a form of a neutral particle 
and scatters taking a form of a charged particle. The 
conversion of a particle is due to interaction with a 
soft background radiation. In the case of an ultrarel- 
ativistic flow, the converter mechanism can be more 
efficient than Fermi acceleration because a neutral 
particle can more easily cross the boundary between 
the flow and the external environment. 

There can be two types of charge exchange cy- 
cles: (1) proton-neutron cycle with proton conversion 
into neutron via photo-meson production and subse- 



quent neutron decay, (2) photon-electron/positron 
cycle with conversion due to photon-photon pair pro- 
duction and inverse Compton scattering. 

The important feature of the second cycle is the 
absence of any limitation of the number of participat- 
ing particle such as the baryon number conservation 
in the first case. High-energy photons and pairs can 
breed exponentially at certain conditions being fed 
by the bulk energy of the flow [8-10]. 

A very simplified description of the photon 
breeding cycle for the case of a relativistic jet can 
be represented by five steps. 

(i) An external high-energy photon of energy e 
enters the jet and interacts with a soft background 
photon producing an electron-positron pair. 

(ii) The time-averaged Lorentz factor of the pro- 
duced pair (as measured in the external frame) gyrat- 
ing in the magnetic field of the jet, becomes 7 ~ T 2 e. 

(iii) The electrons and positrons in the jet Comp- 
tonize soft photons (internal synchrotron or external) 
up to high energies. 

(iv) Some of these photons leave the jet and pro- 
duce pairs in the external environment. 

(v) Pairs gyrate in the magnetic field and Comp- 
tonize soft photons more or less isotropically. Some 
of these Comptonized high-energy photons enter the 
jet again. 

In this cycle, the energy gain, ~ T 2 , is provided 
by the isotropization of the charged particles in the 
jet frame and is taken from the bulk energy of the 
flow. Other steps in the cycle are energy sinks. The 
whole process proceeds in a runaway regime, with 
the total energy in photons and relativistic particles 
increasing exponentially, if the amplification coeffi- 
cient (energy gain in one cycle) is larger than unity. 

This new mechanism is actually much simpler in 
the description than Fermi acceleration, which de- 
pends on detailed geometry of the magnetic field 
and the poorly understood supply of seed nonther- 
mal particles. Once the velocity pattern for the fluid, 
the magnetic field and the external radiation field are 
specified, the fate of each high-energy photon and its 
descendants can be reproduced (in statistical sense) 
from first principles, because the interaction cross- 
section are known with high accuracy. The question 
whether a photon produces a runaway avalanche or 
not, can be answered exactly. 
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Fig. 1. Fraction of the jet power converted into photons ver- 
sus time. The initial density of the high-energy, seed photons 
corresponds to the extragalactic gamma-ray background. Pa- 
rameters: distance from the central engine R = 2xl0 17 cm, jet 
radius Rj = 10 16 cm, jet and disc luminosities Lj = L d = 10 44 
erg s — and Poynting flux of Lb ~ 10 43 erg a -1 . From [10]. 

3. Dissipation of jet bulk energy into 
radiation 

Stern & Poutanen [9, 10] performed a detailed nu- 
merical simulation of the photon breeding in AGN 
jets. The simulation of the particle emission prop- 
agation and interactions was exact and based on 
first principles, while the fluid dynamics was treated 
with a simplified 2D ballistic approximation. Here we 
present some results obtained under assumption that 
the main emission site is the broad emission line re- 
gion [10] . The soft photon background which is neces- 
sary for the converter mechanism was composed from 
direct radiation of accretion disc and an isotropic 
component, which includes reprocessed/scattered ra- 
diation of the disc and surrounding dust. 

To initiate the photon breeding one needs a 
number of seed high-energy photons. Their origin 
is not important as their number can be arbitrarily 
small. Fig. 1 represents a simulation starting from 
the seed gamma-rays corresponding to the extra- 
galactic gamma-ray background. The energy release 
increased by 20 orders of magnitude during 250Rj/c 
(~ 3 years for the jet radius Rj = 10 16 cm) and 
came to the steady state at the level ~ 0.5 of the 
total jet bulk energy. The rapid rise at small times 
corresponds to the exponential growth of the pho- 
ton avalanche as it moves downstream with the flow. 
Then as the avalanche reaches the end of the "simu- 
lation volume" (of length 20Rj) it cannot grow fur- 
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Fig. 2. Working areas of the photon breeding mechanism in 
luminosity-Lorentz factor plane. Solid curve refers to the jet 
power, dashed curves — to the disc luminosity. 

ther, and the growth is further supported by the up- 
streaming photons, which provide a spatial feedback 
loop. 

We have performed several tens simulation runs 
[10] and explored the parameter space, where the 
photon breeding can work. A schematic representa- 
tion of the results is given in Fig. 2. The minimal 
jet Lorentz factor for the supercritical breeding was 
r = 8. At a moderate Lorentz factor the mechanism 
works only at favourable conditions. One of this con- 
ditions is weak magnetic field. The magnetic energy 
flux of the jet should be much smaller than the ac- 
cretion disc luminosity. Otherwise the synchrotron 
losses by pairs at step (iii) of the breeding cycle dom- 
inate over Compton losses. The synchrotron photons 
are too soft to produce pairs and do not participate 
in the breeding, reducing its efficiency. If the jet has 
low magnetization (<r <C 1) or the jet power is smaller 
than the disc luminosity, the breeding is favoured. 

The second favourable condition is the presence 
of a soft isotropic radiation field. The broad line pho- 
tons and scattered photons from the disc are not suf- 
ficiently soft, because pairs in the jet at step (iii) 
interact with them in a deep Klein-Nishina regime. 
Therefore again synchrotron radiation reduces the 
breeding efficiency. A softer radiation can be sup- 
plied by the surrounding dust or by the jet itself. 

The third condition is a high accretion disc lu- 
minosity, Ld > 10 43 /i?i7 erg s _1 (where Rn — 
i?/10 17 cm is the distance from the disc to the site 
of the photon breeding), which is necessary to pro- 
vide sufficient photon-photon opacity across the jet. 
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Fig. 3. Terminal Lorentz factor Tj at the outlet of the cylin- 
dric "simulation volume" with length 20Rj versus distance 
from the jet axis for different initial Lorentz factor of the jet. 

At a higher Lorentz factor the requirement of 
a weak magnetic field is relaxed. If T > 20, then a 
magnetically dominated jet with Lj ~ Ld can radiate 
efficiently. The above conditions are relevant for the 
linear case, when the jet power is moderate and the 
soft synchrotron radiation of the jet is optically thin 
for gamma-rays. The non-linear effects appear at the 
jet power L 3 ; ~ 3 x 10 44 erg s _1 (for R ~ 10 17 cm). At 
Lj > 10 45 erg s _1 the radiative steady-state becomes 
self-supporting even without external radiation. 

The radiative efficiency of the jet increases 
with T and reaches 0.14, 0.56, 0.77, 0.82 for T = 
14,20,30,40, respectively [10]. A high efficiency im- 
plies that the jet undergoes strong deceleration. The 
distribution of terminal T across the jet is shown in 
Fig. 3. One can see that the deceleration is very inho- 
mogeneous: the final Lorentz factor at the jet bound- 
ary is 2-4, while at the center its value is 10-14. 

The resulting spectral energy distributions for a 
sequences of runs with varying luminosity are shown 
in Fig. 4. We see a tendency similar to the observed 
"blazar sequence" : spectra for higher luminosity are 
shifted to lower energy. The observed blazar spec- 
tra demonstrate two distinct components tradition- 
ally interpreted as synchrotron and inverse Compton 
emission peaks of the same electrons. In our sim- 
ulations the synchrotron and Compton components 
are broad and overlap. The possible reason for this 
difference is that the photon breeding mechanism 
works together with other mechanisms like diffusive 
or internal shock acceleration. Moreover, the photon 
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Fig. 4. Simulated spectra for the case Lj = Lj, Lb = 0.2Lj 
and Ld varying from 5 X 10 43 to 10 46 erg s — L 

breeding produces a strong impact on the jet, after 
the rapid inhomogeneous deceleration the jet should 
be highly perturbed. In our interpretation the soft 
blazar component is emitted further downstream due 
to a secondary process, e.g., diffusive reheating of 
cooled pairs produced by the photon breeding mech- 
anism. 
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